Universal decay of scalar turbulence 
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The asymptotic decay of passive scalar fields is solved analytically for the Kraichnan model, where 
the velocity has a short correlation time. At long times, two universality classes are found, both 
characterized by a distribution of the scalar - generally non-Gaussian - with global self-similar 
evolution in time. Analogous behavior is found numerically with a more realistic flow resulting from 
an inverse energy cascade. 
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Recent research on passive scalars in prescribed turbu- 
lent flows has mostly concentrated on universal small- 
scale properties in the presence of large-scale random 
pumping of the scalar (see Ref. |Q and references 
therein). Another important problem is the decay of a 
scalar field 0(r,t) with prescribed random initial condi- 
tions. It is then possible to obtain a stronger form of 
universality at all scales, but only at long times, when 
(most) of the information in the initial conditions has 
been washed out. Specifically, we shall be interested in 
self-similar decay at infinite (or very large) Peclet num- 
bers. We ask the following questions: (i) Is it possible 
to find time-dependent rescaling factors of length and 
scalar amplitude such that in the rescaled variables the 
full single-time probability of the passive field becomes 
independent of time as t — > oo ? (ii) If the first ques- 
tion holds true, what are the universality classes (with 
respect to the initial conditions) of the asymptotically 
self-similar scalar field and what are their statistical prop- 
erties? Such issues are of course not limited to passive 
scalars. For Burgers turbulence, a systematic theory of 
self-similar decay giving, e.g., the law of decay of the 
energy can be derived |2p|| . We shall return to the cor- 
responding issues for 3-D Navier-Stokes turbulence in 
the conclusions. Here we shall first address such issues 
within the framework of the Kraichnan model || of pas- 
sive scalars in which the velocity has a very short corre- 
lation time. This allows analytical and numerical work 
shedding light on the basic mechanisms of scalar decay. 
We show that there are two universality classes, both 
generally having self-similar non-Gaussian distributions. 
Spatially smooth and very rough velocity fields represent 
two limiting special cases: for the former self-similar de- 
cay is broken &i\ while for the latter nearly Gaussian 
self-similar decay obtains. We then turn, using numeri- 
cal simulations, to decay in the presence of more realistic 
velocity fields for a scalar transported by a 2-D velocity 
field arising from an inverse energy cascade process M. 
This also leads to self-similar non-Gaussian decay with 
exponential tails in the single-point scalar pdf. 

The advection-diffusion equation governing the evolu- 
tion of a passive scalar field is 



d t 9{r, t) + v(r, t) ■ V 9(r, t) = nV 2 0(r, t), 



(1) 



where v(r,t) is the incompressible advecting flow and K 
is the molecular diffusivity. Eq. ([!]) governs the evolution 
of the probability density at position r and time t of 
tracer particles. Their Lagrangian trajectories p(t) obey 
the stochastic equation 



dp{t) = v(p(t), t) dt + V2k d/3(t), 



(2) 



where (3(t) is the d-dimensional isotropic Brownian mo- 
tion. The scalar correlation functions at time t can then 
be expressed in terms of the initial ones as 



C N {r,t) 



(6( ri ,t)...6(r N ,t)) 

V N (p,0\r,t) C N (p,0)dp. 



(3) 



Here, r denotes the set r%, . . . , rjv and Vn is the prob- 
ability density that N Lagrangian particles, being at r 
at time t, were previously at p at time t = 0. The ini- 
tial condition 8(r, 0) is taken random and statistically 
homogeneous with correlation length L. The statistical 
homogeneity implies that the absolute position of the 
particles is irrelevant and only the separation variables 
Ti j = Ti — Tj matter. In the sequel we shall always 
consider the propagator in the restricted homogeneous 
sector, i.e. in the (N — l)d separation variables. Note 
that the evolution equation (Q) conserves the so-called 
Corrsin integral || 



J = fc 2 (r,t)dr. 



(4) 



That conservation law alone does not however allow to 
predict the scalar energy decay: for this, we also need 
to know the growth of the scalar correlation length with 
time. This is where the specific properties of the advect- 
ing velocity come into play. 

Let us begin with the Kraichnan model: the velocity is 
statistically homogeneous and stationary, isotropic and 
Gaussian. The second-order correlation function of its 
increments is 

D aP {r, t) = (v a (r, t) Vf}(r, 0)) - (v a (r, t) v p (0, 0)) = 

Dr -i [ r 2 ( 7 + d + 1) 5 aP - (2 - 7 ) r a r p ] S(t) , (5) 
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where 7 controls the roughness of the field. Smooth flows 
correspond to 7 = 0, while for 7 = 2 the opposite limit of 
very rough fields is achieved. Hereafter, we shall concen- 
trate on the generic non-smooth case 7 7^ 0. The growth 
of the scalar length-scale is determined dimensionally to 
be L{t) ~ t 1 / 7 and the units are chosen for convenience 
to set L = D = 1. Note that the velocity field in the 
Kraichnan model is statistically invariant under time in- 
version, i.e. v(r, t) has the same statistical properties 
as —v(r,—t). It follows that the iV-particle propagator 
satisfies 



V N (p, 0\r, t) = V N (p, t\r, 0) = V N (r, t\p, 0) 



(6) 



The major consequence of the short correlation time 
assumption is that the Lagrangian trajectories are a 
Markov stochastic process and that the TV-particle prop- 
agators Vn obey Fokker-Planck equations (see, e.g., 

Ref. My. 



+ -zD a p(r it j) 



V N = 0. (7) 



The advective term in ([!]) conserves all the moments (8 n ). 
The scalar decay is therefore due to the effects of the 
molecular dissipation and, as usual in turbulence, we are 
interested in the singular limit n — > (infinite Peclet 
numbers). The scalar dissipation k((V0) ) remains in- 
deed finite in the limit. From the Lagrangian point of 
view, the molecular noise in ^ can be neglected and the 
singularity of the limit shows up in the fact that initially 
coinciding particles separate in finite time For 
k — ► 0, the fundamental solution of (Q) is self-similar and 
the propagator possesses the scaling property 

V N (Xp,0\Xr,X^t) = X^ N ^ d V N (p,0\r,t) . (8) 

The simplest objects to investigate are the single-point 
moments (9 2n )(t) and we are interested in their long-time 
behavior t > 1. The first universality class is formed by 
the initial conditions where the Corrsin invariant Jo de- 
fined in (^) does not vanish. We shall also make the nat- 
ural assumption that the initial scalar correlations decay 
rapidly enough to ensure integrability of the cumulants 
of arbitrary order. The long-time behavior of the single- 
point moments is obtained by setting r = in (^) and 
using (||) to rescale out the time dependence. The re- 
sulting object involves t d ' 7 C2(pt 1 ' 7 , 0) which tends to 
S(p)Jo- The final result is 

jn 

(^ 2 ")(*)- ^7(2^-1)'! 

x J V 2n (Pi, Pi, ■ ■ ■ Pn 1 Pn 1 

1|0,0) d Pl ...dp n , (9) 

as is easily seen for the special case of initial Gaussian 
statistics. In fact, it is not hard to show that the same 




FIG. 1. The normalized kurtosis (6< 4 )/3{6> 2 ) 2 vs time in the 
Kraichnan model for the two values 7 = 5/4 (dashed line) 
and 7 = 2/3 (solid line). The error bars are comparable to 
the size of the symbols. 

result holds for all initial data in the universality class 
Jo 7^ and rapid correlation decay. The integrability 
condition on the cumulant of order 2n implies indeed 
that it must depend on all the (2n — 1) independent in- 
terparticle separations. In the limit of long times, each 
dependence brings a function S(p i j) and a factor t d ^ in 
the denominator. The total non-Gaussian contribution 
from the initial conditions will therefore decay faster than 
t~ nd h . Note that this argument implies a stronger form 
of universality than in forced steady states: the asymp- 
totic pdf is here independent of the initial data, whereas 
in the steady states the scalar pdf depends on the details 
of the force. The expression (^) shows that the scalar 
statistics is self-similar in time. For example, the sin- 
gle point pdf V(6,t) takes the form t d / 2 ^Q{9t d / 2 ~<). The 
propagator in the limit of initially coinciding points that 
appears in (^) is a regular function of its arguments p0| . 
The smooth case 7 = is special in this respect since 
nearby particles separate exponentially and their separa- 
tion is proportional to its initial value. This is the phys- 
ical reason why @ and the ensuing self-similarity of the 
scalar distribution do not hold for 7 = (see Refs. |^,0|). 
Note that the asymptotic pdf V(8,t), although univer- 
sal and self-similar, is generally non-Gaussian. From (|]), 
the hypothesis of Gaussianity is equivalent to the factor- 
ization of the probability for the In particles to collapse 
into pairs, i.e. the integral in (|^) should be equal to 
[/ V2 (p, P, 1 |Q, 0) dp] . This is generally not the case, 
due to the correlations existing among the particle tra- 
jectories. The only exception is 7 = 2 where the par- 
ticles are independent. The degree of non-Gaussianity 
is thus expected to decrease when 7 increases, as con- 
firmed in Fig. 1. The simulations were performed by the 
Lagrangian method presented in Refs. pj]Jl^ |. 

The second universality class is for initial conditions 
with Corrsin invariant Jo = 0. From the realizability of 
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FIG. 2. The scalar energy (6 2 ) vs t for the advection by a 
2D flow generated by an inverse energy cascade. The dashed 
line is the asymptotic power law behavior 1/t 3 . 

the field it follows that the integral J\ = jr' y C2(r)dr 
is non- vanishing jljj. It is also easily checked using (||) 
and (0) that the integral remains invariant in time. Jo = 
implies that the integral of the correlation CV with 
respect to any of its arguments vanishes. Furthermore, 
the single-point moments can then be expressed as 



l )(t) = J Q 2 „ (p,t\r,0) C 2n (p,0) dp. 



(10) 



Here 



2n 



Q2n (p,t\r,0) = \J[ A p j S 2n V 2n (p,t\r,0) , (11) 

S2n denotes the symmetrizer in the p variables and 
the difference operator A acts on a function f(x) as 
A x ,yf(x) =f(x)-f(y). To obtain @ , we set r = in 
(||), symmetrize in the p variables and use the condition 
Jo = to reexpress the result in terms of difference oper- 
ators. For this, note that whenever any of the p variables 
in the propagator is replaced by the integral vanishes. 
The long-time behavior of the single-point moments is 
then derived from ( |Io|) exploiting again the scaling prop- 
erty m). For Gaussian initial conditions, we obtain 



jn 

(OW-^(^ T Ty(2n-l)!!x 

0.2-a (Pi, P 2 , ■ ■ ■ Pln-liPlm l|fi, °) 



(12) 



lim 



[P 



1,2 ■ 



• P2n-l,2nJ 



H d P2 



i=l 



where we used the convergence of A d / 7 (Ap) 7 C 2 (Ap, 0) to 
8(p)J\ for large A's. The limit of ■ ■ ■ P2n-i.2n tending 
to zero, which is understood in (p"2|), has been shown to 
be well defined |l4|] . The strong universality of the scalar 
distribution is proved as for the Jo ^ case. The ex- 
pression ( |l2| ) indicates that the asymptotic scalar distri- 
bution is again self-similar although the rescaling factor 



is now i n agreement with the exact solution for 

the second-order case Q. 

The long-time behavior of multi-point scalar correla- 
tions can be derived from (^) as for the single-point mo- 
ments. Using the scaling property (||) we are generally 
led to analyze Vn{p, lln/i 1 ^ 7 , 0), i.e. the limit where the 
N particles are initially close. This limit is crucial also 
for the small-scale scaling properties in the presence of 
pumping and leads to an asymptotic expansion for the 
propagator [jl(| 

V N (P,t\*Z,0) ~ Yl XPa+7k 5a,fe(pJ), (13) 



valid for small A. The first sum is over the zero modes 
f a = f a0 with scaling dimensions p a (in increasing or- 
der). The functions f a ^ for different k form the corre- 
sponding towers of the so-called slow modes. The lead- 
ing term in the expansion ( |lg| ) comes from the constant 
zero mode /o.o = 1 which corresponds to gofi(p 1 t) — 
'Pn (p Jlfij 0), i.e. the propagator for N particles that 
are initially at the same position. Inserting ( p"3| ) into 
(0) and using the dominance of go,o is is easy to check 
that generic multi-point correlations will asymptotically 
decay as the corresponding single-point moment. Sub- 
dominant terms emerge by combining the correlation 
functions to cancel the dominant contribution. The 
classical example is given by the structure functions 
S 2n (r,t) = ([6(r,t) - 6»(0,t)] 2n ), related to the corre- 
lation functions by 



S 2n {r,t) 



/ 2n \ 

Y[A ri0 \C 2n (r,t) 



(14) 



The difference operator A was defined in (|ll|). The cru- 
cial point to remark is that all the reducible contributions 
to the correlations, i.e. not depending on all the coordi- 
nates n, . . . r2n, are annihilated by the combination of 
the structure functions. Inserting ( |l3| ) into (||) we obtain 



S 2n 



(n^lA r4 , )/2n,0(r.) 



P2n 



52n,o(£, 1) C 2 „(^p,0) dp 
(15) 



where /2n,o is the irreducible zero mode in ([13|) with the 
lowest dimension. The quantities Q 2n — p 2 n coincide with 
the exponents of the structure functions in forced steady 
states |l(J. This justifies the classical notion of quasi- 
cquilibrium, i.e. that short-distance scaling laws are the 
same in forced and decaying situations provided the basic 
parameters are made appropriately time-dependent. 

The analytical solution found for the Kraichnan model 
permits to capture the basic features of the scalar de- 
cay. It is in particular quite clear that the key ingredient 
for the asymptotic self-similarity of the scalar pdf is the 
validity of the rescaling property m). We also performed 
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mittent spatial distribution of the kinetic energy? It is of 
great interest to address such issues experimentally and 
by advanced direct numerical simulations. 
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FIG. 3. The rescaled scalar pdf's vs t z 6 for the inverse en- 
ergy cascade flow at the three times marked by a star in Fig. 2. 
The values on the abscissa are normalized by their r.m.s. val- 
ues. 

simulations using a more realistic "inverse energy cas- 
cade" 2-D velocity field which show that such property 
is not limited to short-correlated velocities and is ex- 
pected to hold (at least for large enough times) for any 
self-similar velocity field. Indeed the "inverse cascade" 
flow is scale- invariant (no intermittency) of exponent 1/3, 
isotropic and, most important, has finite realistic corre- 
lation times [jl5|— p~T|] . The scalar initial distribution is 
taken Gaussian with correlations decaying exponentially 
and a correlation length comparable to the velocity in- 
jection scale (see Ref. jl7j for more details). In Fig. 2 we 
present the curve of the energy decay showing that the 
exponent coincides with the value d/j obtained from the 
Kraichnan model (7 = 2/3 for the inverse cascade flow). 
The rescaled scalar pdf's V(0, t)/t 3 for three times in the 
asymptotic power-law region are presented in Fig. 3. It 
is quite clear that the asymptotic distribution is again 
self-similar and its tails are exponential. The tails are in- 
deed dominated by the rare events governed by Poisson 
statistics where the scalar is weakly mixed by diffusive 
effects @. 

For passive scalar decay in an intermittent velocity 
field, such as a 3-D turbulent flow, the presence of various 
scaling exponents makes it unlikely that the propagator 
possesses a rescaling property like (|8|). The self-similarity 
of the scalar decay found here might then be destroyed, 
with the kurtosis and higher-order scalar moments pos- 
sibly growing indefinitely with time. An even more chal- 
lenging problem is that of the decay of 3-D turbulent 
flow which has been frequently addressed, but usually in 
a phenomenological way and with the emphasis on the 
law of energy decay (see Refs. for a review). Is 

there self-similar and universal decay as for the problem 
considered in this paper or is the kurtosis of the velocity 
growing indefinitely, possibly leading to a highly inter- 
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